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ABSTRACT 

We report a series of high-resolution cosmological N-body simulations designed to explore the 
formation and properties of dark matter halos with masses close to the damping scale of the primordial 
power spectrum of density fluctuations. We further investigate the effect that the addition of a random 
component, Wrms, hito the particle velocity field has on the structure of halos. We adopted as a fiducial 
model the A warm dark matter (AWDM) cosmology with a non-thermal sterile neutrino mass of 0.5 
keV. The filtering mass corresponds then to Mf = 2.6 x IO^^/i^^Mq. Halos of masses close to Mf 
were simulated with several million of particles. The results show that, on one hand, the inner density 
slope of these halos (at radii ^ 0.02 the virial radius R„) is systematically steeper than the one 
corresponding to the Navarro- Frenk- White (NFW) fit or to the cold dark matter counterpart. On the 
other hand, the overall density profile (radii larger than 0.02 R^,) is less curved and less concentrated 
than the NFW fit, with an outer slope shallower than -3. For simulations with Wims, the inner halo 
density profiles flatten significantly at radii smaller than 2-3 /i~^kpc (^ O.OlORi, — O.OlSRi,). A 
constant density core is not detected in our simulations, with the exception of one halo for which the 
flat core radius is w l/i~^kpc. Nevertheless, if "cored" density profiles are used to fit the halo profiles, 
the inferred core radii are k, (0.1 — 0.8)ft.~^kpc, in rough agreement with theoretical predictions based 
on phase-space constrains, and on dynamical models of warm gravitational collapse. A reduction of 
frms by a factor of 3 produces a modest decrease in core radii, less than a factor of 1.5. We discuss the 
extension of our results into several contexts, for example, to the structure of the cold DM micro-halos 
at the damping scale of this model. 

Subject headings: cosmology: dark matter — galaxies:halos — methods:N-body simulations 



1. INTRODUCTION 

The nature of dark matter (DM) is one of the most 
intriguing and fundamental problems in cosmology and 
particle physics. The standard hypothesis assumes that 
dark matter is made of non-baryonic collisionless ele- 
mental particles that become non-relativistic very early 
in the history of the Universe (cold). This minimal sce- 
nario, named Cold Dark Matter (CDM), has successfully 
explained the observed structure of the universe at large 
scales, like the two-point correlation function of galax- 
ies and the Cosmic Microwave Background Radiation 
(CMBR) anisotropics (see for recent results Springel et 
al. 2006; Spergel et al. 2006). Confrontation of model 
predictions with obervations turns out to be more com- 
plicated at galactic scales, because non-linear dynamics 
and baryonic processes may distort considerably the un- 
derlaying DM distribution. Thus, the predictions of the 
CDM scenario and its variants at the scale of galaxies, 
are an active subject of study. 

Two of the most controversial CDM predictions are: 
(i) the large abundance of subhalos in galaxy-sized halos 
(Kauffmann et al. 1993), and (ii) the cuspy inner density 
profile of dark halos (e.g., Navarro et al. 2004; Diemand 
et al. 2004). Based on comparison with observations, it 
has been suggested that both predictions may indicate a 
flaw of the CDM scenario (Klypin et al. 1999a; Moore 



et al. 1999a; Moore 1994; de Blok et al. 2001, Gentile et 
al. 2005,2007). However, these comparisons might be bi- 
ased by astrophysical processes in action during galaxy 
assembly and evolution, as for example the inhibition 
of star formation in small subhalos (Bullock, Kravtsov, 
& Weinberg 2000; Benson et al. 2003; Governato et al. 
2007) or the halo core expansion due to energy or angular 
momentum transfer from dark/baryonic structures (Ma 
& Boylan-Kolchin 2004; El-Zant et al. 2004; Weinberg 
& Katz 2007), but see also (Colin et al. 2006; Ceverino 
& Klypin 2007; Sellwood 2007). It has also been showed 
that the disagreements may be a consequence of system- 
atics in the observational inferences, (e.g., Simon & Geha 
2007; Rhee et al. 2004; Hayashi & Navarro 2006; Valen- 
zuela et al. 2007). 

It is also possible that slight modifications to the CDM 
particle properties could solve or ameloriate the men- 
tioned potential problems if they persist (e.g., Spergel 
& Steinhardt 2000; Colin et al. 2002). As the precision 
of observations and the control of systematics improve, 
the confrontation with model predictions opens a valu- 
able window for constraining the dark matter properties. 
Among the "slight" modifications with the CDM sce- 
nario is the introduction of warm particles, instead of 
cold ones (Warm Dark Matter, hereafter WDM). WDM 
implies two more degrees of freedom in comparison to 
CDM: (i) a damping (filtering) of the power spectrum at 
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some intermediate seale, Aj, due to the rclativistic free 
streaming, and (ii) some primordial random velocity in 
the dark particles, Vrms (for CDM, the mass correspond- 
ing to A/, Mf, is comparable to planet masses, e.g., Die- 
mand et al. 2005; Profumo et al. 2006, and Wrms ~ 0). 

Previous numerical studies have shown that WDM can 
be very effective in reducing the amount of (sub)structure 
below the filtering mass Mf (e.g., Colfn, Avila- Reese, & 
Valenzuela 2000; Avila-Reese et al. 200f; Bode et al. 
2001; Knebe et al. 2002). In addition, both the peculiar 
dynamical formation history of these halos and its ran- 
dom Wrnis can also impose an upper limit on the phase 
space density, potentially producing an observable core 
of constant density (Hogan & Dalcanton 2000; Avila- 
Reese et al. 2001). If Vrms bas a thermal origin its am- 
plitude is linked directly to the mass of the WDM par- 
ticle (Hogan & Dalcanton 2000); notice, however, that 
the amplitude of the random velocities may depend on 
other physical factors not directly related to the parti- 
cle mass. This is the case, for instance, of gravitinos 
produced non-thermally by late decays of the Next to 
Lighest Supersymmetric Particle (NLSP) (Feng, Rajara- 
man & Takayama 2003; Strigari et al. 2007, see for a 
recent review Steffen 2006). Thus, an exploration of the 
effect of a random velocity component independently of 
the dark matter particle mass seems to be necessary. 

Predictions of the core radius in WDM halos have been 
computed assuming a King profile (Hogan & Dalcanton 
2000) or a subclass of the Zhao (1996) profiles (Strigari 
et al. 2006, herafter S2006). It is still unknown which 
estimates give the more accurate value, yet these pre- 
dictions are necessary for comparison with observations 
(see e.g., S2006; Gilmore et al. 2007). Unfortunately, 
the predicted core radii for masses of the most popular 
WDM particle, the sterile neutrino, allowed by observa- 
tional constraints (> 2 keV, see §5 for references) are 
below the resolved scales in current simulations. 

The structure of WDM halos of scales below the filter- 
ing radius A / might be different from their CDM coun- 
terparts not only in the central parts but also in their 
overall mass distribution. This is somehow expected be- 
cause the assembly history of these halos is different from 
the hierarchical one. Besides, they form later and have 
concentrations lower than those ones derived for CDM 
halos (Avila-Reese et al. 2001; Bode et al. 2001; Knebe et 
al. 2002). However, it is controversial whether the shape 
of the density profile differs systematically from the cor- 
responding CDM one (see for different results e.g., Huss 
et al. 1999; Moore et al. 1999b; Knebe et al. 2003). On 
the other hand, WDM halos of mass close to Mf (with 
Vrms set to 0) Can be thought of as scaled-up versions of 
the first microhalos in a CDM cosmology. The earliest 
collapse of CDM microhalos is a subject of considerable 
current interest (e.g., Diemand et al. 2005; Gao et al. 
2005). 

In this paper we explore by means of numerical sim- 
ulations the two questions mentioned above: the overall 
structure of dark halos with masses close or just below 
the cutoff mass in the power spectrum of fluctuations, 
and the inner density profile of these halos when a ran- 
dom velocity is added to the particles. For our numerical 
study we use the truncated power spectrum correspond- 
ing to a non-thermal sterile neutrino of mw = 0.5 keV 
{Mf = 2.6 X 10^2 /i-IMq). We initially neglect the parti- 



cle random velocity (wrms=0), and later we will consider 
two values of Vj-^g that cover the range of velocities cor- 
responding to thermal and non-thermal mw = 0.5 keV 
WDM particles. It is important to remark that our goal 
is not to study a specific WDM model but to explore 
in general the effects on halos of the power spectrum 
truncation and the addition of random velocities to the 
particles. 

The structure of the paper is as follows. In §2 we de- 
scribe the cosmological model that we use for our inves- 
tigation : a WDM model with a filtering radius at the 
scale of Milky Way size halos. Halos of these scales were 
resimulated with higher resolutions, first without adding 
the corresponding random velocity, and then with adding 
this velocity component to the particles. Details of the 
numerical simulations carried out in this paper are given 
in §3. The results from our different simulations are pre- 
sented in §4. Section §5 is devoted to a discussion of the 
results and their implications. Finally, in §6 we present 
the main conclusions of the paper. 

2. THE COSMOLOGICAL MODEL 

The general cosmological background that we use for 
our numerical simulations corresponds to the popular flat 
low -density model with Hq = 0.3, Ua = 0.7 and h = 0.7 
(the Hubble constant in units of 100 km sec^^ Mpc~^). 

For the experiments designed to explore the struc- 
ture of dark matter halos with masses close to Mf, we 
adopt an initial power spectrum corresponding to a non- 
thermal sterile neutrino of 0.5 keV. Even if this WDM 
model is ruled out by observations (see §5 for references), 
it is however adequate for the purposes stated in the In- 
troduction. As is shown below, the filtering mass Mf 
corresponding to this WDM particle is of the order of 
Milky Way-size halos, namely the halos that we are able 
to follow with high resolution in a cosmological simula- 
tion. The high resolution of the simulations avoids that 
the formation of halos with a mass scale near to Mf will 
be dominated by discreteness effects (see §3.1). Besides, 
for the resolution that we attain, we expect to resolve the 
inner regions where a flattening in the halo is predicted 
density profile is predicted for the case when Urmsis in- 
troduced. 

We use here the transfer function Tg for the non- 
thermal sterile neutrino derived in Abazajian (2006a). 
The WDM power spectrum is then given by 

PwDM{k)=T^{k)PcDM{k), (1) 



where 



T,ik) = [1 + {akry 



(2) 



and PcDM is the CDM power spectrum given by Klypin 
& Holtzman (1997). This fit is in excellent agreement, 

at the scales of interest, with the power spectrum ob- 
tained with linger, which is contained inside the cosmics 
package^. The parameter a is related to the mass of the 
sterile neutrino, 0dm and h through 



DM 



VlkeVy V 0.26 



n 



-] /.-Mpc, 



(3) 



where a = 0.189, b = -0.858, c = -0.136, d = 0.692, v = 
2.25, and /x = 3.08. The power spectrum is normalized 

^ http: / / web.mit.edu/edbert /cosmics- 1 .04.taj.gz 
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to CTg ~ 0.8, a value close to that estimated from the 
third- release of WMAP (Spcrgcl et al. 2006). Here cts 
is the rms of mass fluctations estimated with the top-hat 
window of radius 8/i~^Mpc. 

As in Avila-Reese et al. (2001), we have defined the 
free-streaming (damping) wavenumber, fc/, as the k for 
which the WDM transfer function r^(fc) decreased to 
0.5, and compute the corresponding filtering mass in the 
linear power spectrum as 

where p is the present day moan density of the universe. 
The filtering wavelength A/ = li^jk^ of a non thermal 
sterile neutrino of 0.5 keV is 3.9 /!.~^Mpc, which corre- 
sponds to a filtering mass Mj = 2.6 x 10^^ /i^^Mq. The 
random component was linearly added to the velocities 
calculated with the Zel'dovich approximation at the on- 
set of the simulation. 

A number of Milky Way-size halos are simulated with 
the WDM power spectrum given by eq. (1) and with 
Vrms=0- In this way we isolate the effects of the power 
spectrum filtering on the structure of the halos. Later, 
the same halos are resimulated but adding random veloc- 
ities to the particles. We approximate the shape of the 
particle phase space distribution function (DF) with the 
corresponding thermal one; i.e, we use a thermal equi- 
librium Fcrmi-Dirac DF. This is a good approximation 
for a non thermal sterile neutrino of 0.5 keV (Abazajian 
2006a). Further, as a first approximation, we assume 
that the amplitude of the DF of our non-thermal ster- 
ile neutrinos is equal to the one corresponding to their 
thermal partners with mw=0.5 keV (Hogan & Dalcanton 
2000). With this assumption for DF shape and ampli- 
tude, the neutrino random velocity at z = 40, the initial 
redshift in our simulations, is « 4 km/s. 

The particle velocity dispersion in the linear regime 
decreases adiabatically with the expansion. The conven- 
tion is to define Wrms in terms of the value extrapolated 
to the present epoch, 'yrms(z) = Vrms(0)(l-l--2:); therefore, 
for our case frms(O) « 0.1 km/s. We should note that 
the rms velocity amplitudes of sterile neutrinos and ther- 
mal particles of the same mass can not be the same, but 
larger because the former become non-relativistic later 
than the latter. A rough calculation shows that for the 
0.5 keV sterile neutrino, Vims should be approximately 
three times higher than for a 0.5 keV thermal WDM par- 
ticle, i.e. Wrms(O) ~ 0.3 km/s. We resimulate halos with 
both values of the random velocity, Wrms(O) ~ 0.1 and 0.3 
km/s, having in mind that our main goal is not to study 
a particular WDM model but to explore the structure of 
the dark halos with different initial conditions . 

3. NUMERICAL SIMULATIONS 

A series of high resolution simulations of Milky Way 
size galactic halos are performed using the Adaptive Re- 
finement Tree (ART) N-body code (Kravtsov, Klypin, & 
Khokhlov 1997) in its multiple mass scheme (Klypin et 
al. 2001). The ART code achieves high spatial resolu- 
tion by refining the root grid in high-density regions with 
an automated refinement algorithm. 

In all of our experiments we first run a low-mass resolu- 
tion (LMR) simulation of a box of L^ox = 10/i~^Mpc on 



a side and 64'^ or 128"^ particles. We then select a spher- 
ical region centered on a Milky Way size halo of radius 
~ 3 times the virial radius^ of the chosen halo. The 
Lagrangian region corresponding to the z = spherical 
volume is identified at 2; = 40 and resampled with ad- 
ditional small-scale waves (Klypin et al. 2001). Halo 
names (first column in Table 1) are denoted by a capital 
letter followed by the effective number of particles that 
were used in the high resolution zone and a subindex 
that indicates the value of Wrms(O) (in km/s; reported in 
column 2)^. The mass per particle, nip, in the high- 
resolution region is given in column (5) while in column 
(6) we give the virial mass of the halo, M^,. The latter 
along with rrip can be used to compute the number of 
particles inside a halo. The global expansion timestep 
Aao, and the formal spatial resolution hfor -measured 
by the size of a coll in the finest refinement grid arc 
given in columns (3) and (4). ART integrates the equa- 
tions using the expansion factor a as the time variable 
such that at 2; = 0, a = 1. 

We started the simulations at 2; = 40 because the power 
at the Nyquist frequency at this redshift is in the linear 
regime. Note that spurious noise can influence the evo- 
lution of simulations that include Wrms if they start too 
early. The noise might be particularly important when 
the generated ?;rms have amplitudes comparable to those 
of the Zel'dovich peculiar velocities. In order to show 
this, we started a simulation aX z = 100 and evolved it 
up to z = 40. At 2 = 100 the v-cms velocities are on av- 
erage 2.5 times greater than at 2: = 40. Figure 1 shows 
the measured initial power spectrum of the simulation 
started at z = 40 (solid line) and the power spectrum 
for the simulation started at z = 100 and measured at 
2; = 40 (squares). The latter simulation developed spuri- 
ous power at frequencies higher than about 0.6 /i/Mpc. 
We run another simulation started at z = 100, but with 
noaddition of a relic velocity component. In this case 
we did not detect the evolution of spurious power by 
z = 40. In order to make sure that this spurious noise 
does not appear when the simulation is started at z = 40 
we repeated the experiment but now the onset of the 
simulation is set at z = 40 and the power spectra are 
measured at z = 20. Unlike the previous case, no differ- 
c!ncc!s bctwcicn the power spectra were detected. In other 
words, as far as the initial power spectrum is concerned 
it does not matter if the simulation is started at z = 40 
o z = 20. 

Concern may arise about the structure of halos simu- 
lated in a relatively small Lbo,i; = 10/?,^^Mpc box, spe- 
cially in a WDM cosmology, where there is a scale below 
which the power spectrum exponentially drops to zero. 
Avila-Reese et al. (2001) discussed this potential issue 
and concluded that in order to be confident about the 
simulated halo structure, a box size greater than the fil- 
tering length A/ should be used. For our simulations 
Lftox is 2.5 larger than A/. In any case, we also exper- 
imented with other box sizes, namely, ISft.^'^Mpc and 
20ft.~'^Mpc (not shown in the Tables), for the v^^s = 

^ This radius is defined as the radius at which the mean overden- 
sity is equal to the predicted by a top-hat spherical collapse which 
is 337 for our selected cosmological model. 

The random velocity component was linearly added to the 
velocities correspondingt to the Zeldovich approximation at the 
onset of the simulation. 
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TABLE 1 
Model and Halo Parameters" 



name tag 


Urms 


Aao 


hfor 


nip 


Mv 


Knax 


R.7J 


Cl/5 


cnfw 


Lgpi% 


rc 




c,30 




km/s 


10-3 


kpc/h 


lO^^Mo/fc 




km/s 


kpc/h 






f/i^Mm/pc^l 

L KZJ I y J 


kpc/h 


kpc/h 


kpc /h 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


A2560.0 


0.0 


1.0 


0.305 


4.96 


2.84 


224.1 


286 


5.6 


8.3 


-1.44 


- 


- 


- 


A256o.i 


0.1 


1.0 


0.305 


4.96 


2.37 


218.2 


228 


5.0 


5.9 


-1.44 


0.082 


1.151 


0.572 


B2560.0 


0.0 


0.5 


0.305 


4.96 


4.36 


253.2 


330 


5.1 


7.7 


-1.46 








B256o.i 


0.1 


0.5 


0.305 


4.96 


4.35 


250.8 


330 


5.0 


6.7 


-1.54 


0.098 


1.380 


0.565 




U.U 


U.o 


yj. i-oZ 


yj.DZ 


i.zy 


1 fin 7 


zzz 




fl A 
D.4 










C512o.i 


0.1 


0.5 


0.152 


0.62 


1.26 


158.9 


220 


4.3 


4.7'' 


-1.79 


0.093 


1.265 


0.689 


C5120.3 


0.3 


0.5 


0.152 


0.62 


1.28 


158.5 


221 


4.3 


3.6'' 


-1.99 


0.121 


1.527 


0.827 


D5120.0 


0.0 


0.5 


0.152 


0.62 


1.28 


160.1 


221 


4.3 


6.1 


-1.54 








D5120.1 


0.1 


0.5 


0.152 


0.62 


1.24 


160.9 


219 


4.6 


4.2*' 


-1.89 


0.106 


1.389 


0.816 


D512o lb 


0.1 


0.5 


0.152 


0.62 


1.25 


159.6 


219 


4.3 


4.8'' 


-1.77 


0.091 


1.246 


0.598 


D512o.3b 


0.3 


0.5 


0.152 


0.62 


1.26 


160.4 


220 


4.3 


3.7 


-1.96 


0.118 


1.498 


0.792 


E5120.0 


0.0 


1.0 


0.152 


0.62 


2.13 


197.3 


262 


4.9 


6.4 


-1.57 








E512o.i 


0.1 


1.0 


0.152 


0.62 


2.06 


192.3 


259 


4.7 


5.7*' 


-1.66 


0.092 


1.287 


0.657 


E5120.3 


0.3 


1.0 


0.152 


0.62 


2.14 


197.0 


262 


4.9 


5.3'' 


-1.71 


0.101 


1.380 


0.519 



^AU the halos presented in this Table were resimulated from a L^o^ lO/i ^Mpc box simulation. 

''Note that the NEW function does not provide a good fit to the density profiles of halos with frms > 0. However, for completeness, 
we report here the value of cnfw obtained from the fit. 




0.5 1 1.5 



log k(h/Mpc) 

ElG. 1. — Comparison of the power spectra measured at z = 40 
for the simulation started at z = 40 (solid line) and the one started 
at 2 = 100 (squares). The longest plotted wavelength is L(,oa; while 
the highest frequency is {2it / Lf^^x) 256. 

case, and found results similar to those reported here for 
the lO/i-^Mpc box. 

The bound density maxima (BDM) group finding algo- 
rithm (Klypin et al. 1999b), or a variant of it (Kravtsov 
et al. 2004) , is used to locate the halos in the simulations 
and to generate their density profiles. The BDM finds 
positions of local maxima in the density field smoothed 
at the scale of interest and applies physically motivated 
criteria to test whether a group of particles is a gravita- 
tionally bound halo. 

Aside from those halos shown in Table 1, for the halo 
D with Wrms=0 we have also run a very high resolution 
simulation with about 31 million particles in the high 



resolution zone. This halo was taken from the same 10 
ft.~^Mpc box run but with 128'^ particles in the LMR 
mode. As far as we know this is the highest resolution 
simulation of a halo run in a WDM cosmology. The same 
halo was also simulated with less resolution for a conver- 
gence test. The parameters of the sequence of halos D 
are resumed in Table 2. 

The simulations presented here differ in several aspects 
from previous WDM simulations. First, it should be em- 
phasized that our aim rather than discussing a specific 
WDM model is to explore the influence of the truncation 
of the power spectrum and/or the addition of random ve- 
locities upon the structure of dark halos of masses close 
to the truncation scale. For this aim we need to simu- 
late (a) halos with very high resolution, and (b) halos 
with masses close to Mf. The halos simulated in Avila- 
Reese et al. (2001) had several times less particles than 
the best-resolved halos presented here and the aims in 
that paper focused in exploring general halo properties 
for a concrete WDM model. Other papers aimed to study 
the properties of WDM halos (Bode et al. 2001; Knebe 
et al. 2002; Busha et al. 2006) focused more in the sta- 
tistical aspects than in details of the inner halo struc- 
ture; therefore, the halos in these papers had resolutions 
much lower than those attained here. The properties of 
the WDM halos simulated here are in general agreement 
with previous findings; for example, their concentrations 
are systematically lower (Avila- Reese et al. 2001; Eke et 
al. 2001; Bode et al. 2001) and they form later (Knebe 
et al. 2002; Busha et al. 2006) than the corresponding 
ACDM halos. 

3.1. Discreteness effects 

One of the motivations of this paper is to investigate 
the structure of well resolved halos with masses close or 
below the damping (truncation) scale in the power spec- 
trum, Mf. The origin of these halos is controversial. 
Halos with masses close to Mj (truncation halos) could 
be formed by a quasi-monolithic collapse of filaments of 
size ~ A/ (e.g., Avila-Reese et al. 2001). They could also 
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TABLE 2 

High Resolution Halo D Parameters 



^hox 


name tag 


^rms 


timcstcp 


resolution 


nip 




(ft-iMpc) 






(10-3) 




(A-IM0) 


(IOI2 /i-1Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


10 


D10240.0 


oflf 


0.5 


0.040 


7.75 X 10* 


1.28 


10 


D5120.0 


oflf 


0.5 


0.152 


6.20 X 10^ 


1.27 


10 


D2560.0 


oflf 


0.5 


0.305 


4.96 X l(f 


1.27 


10 


D1280.0 


oflf 


0.5 


0.610 


3.97 X 10^ 


1.25 



be just the result of an incomplete collapse, highly de- 
viated from the spherical-symmetric case, of originally 
larger structures assembled hierarchically (Busha et al. 
2006). On the other hand, it has been suggested that ha- 
los with masses considerably less than M f form by frag- 
mentation of the shrinking filaments of size ~ A/ (e.g., 
VaUnia et al. 1997; Avila- Reese et al. 2001; Bode et al. 
2001; Gotz & Sommer-Larsen 2003; Knebe et al. 2003). 
However, it is also known that the filaments in Hot Dark 
Matter simulations that start from a cubic lattice break 
up into regularly spaced clumps, which reflect the ini- 
tial grid pattern. Therefore, some of these halos seen 
in WDM simulations could be spurious, product of dis- 
creteness effects. Recently, Wang & White (2007) have 
shown that this artifact is present even for a glass-like 
initial particle load (White 1996). 

As Wang & White (2007) show, halos of masses smaller 
than a given effective fraction of M/, which depends on 
the resolution of the simulation, will be spurious. We se- 
lected the WDM model (§2) and the number of particles 
in the simulations in such a way that the halos studied 
here can neiher be fake nor affected by discreteness ef- 
fects. Following Wang & White (2007), in our model 
with 1024^ effective number of particles, only structures 
with masses lower than 1.3 xlO'^/i-^M© arc candidate to 
be spurious. For the models with 512"^ and 128'^ effective 
number of particles, the masses of structures triggered 
by the initial grid spacing are < 2.6 x 10^h~^MQ and 
1.04 xlO^°/i~^M0, respectively. For our WDM model, 
the first halos to form have masses > lO^^ft^^M©, which 
are far from the mentioned effective resolution limits. 

4. RESULTS 

We first present the results of our WDM simulations 

with Wring = 0. In this way, we explore the inner structure 
of halos close to the truncation mass Mf, which could be 
scaled-up versions of the first CDM microhalos (of masses 
10~^Mq for a neutralino mass of 100 GeV). After- 
wards we present the results of the same simulations but 
introducing random velocities to the particles with two 
different amplitudes, 'yrms(O) = 0.1 and 0.3 km/s (see 
§2). The main goal of the last simulations is to explore 
the predicted flattening in the halo inner density profile 
produced by the addition of a random component to the 
particle velocities(see section 1). Table 1 resumes the 
main properties of the resimulated halos, which were se- 
lected to have masses close to the truncation mass of the 
initial power spectrum. 

4.1. The Structure of Halos at the scale of Damping 

Figures 2 and 3 show the spherically averaged den- 
sity profiles measured for the halos of masses around the 



power-spectrum filtering mass, Mf, and with t;ruis=0. 
The first panel of Fig. 2 shows halos A256o.o and 
B256o.o (the latter was shifted by —1 in the log): the sec- 
ond and third panels show halos C512o.o and E512o.o, re- 
spectively. In Fig. 3, the same halo D but simulated with 
4 different resolutions is shown. In both Figs., the thin 
dashed lines are the best Navarro- White-Prenk (NFW, 
Navarro, Frenk, & White 1997) fits to the showed den- 
sity profiles. In the lower panels, the residuals of the 
measured density profile and the NFW fit are plotted 
with the same line coding as in the corresponding upper 
panels. 

Halos A256o.o and B256o.o in the left panel and 
C512o.o and E512o.o in the second and third panels 
of Fig. 2, were simulated with formal spatial resolutions 
hfor = 0.305/i^^kpc and 0.152/i-^kpc, respectively (see 
Table 1). Previous convergence studies for CDM halos 
simulated with the ART code have shown that the in- 
nermost halo density is reliable only for radii larger than 
four times /ifor and containing more than 200 particles 
(Klypin et al. 2001). For all the density profiles shown 
in Figs. 2 and 3, the innermost plotted point corresponds 
to radii larger than /i4 = 4 x /ifor by ~ 30% and they con- 
tain more than 200 particles. The convergence analysis 
that we have carried out for our WDM halos suggests 
that instead of /14, the innermost radius should be close 

to 8 X hfor, hs- 

Figure 3 compares the density profiles of halo D, which 
was re-simulated with four different resolutions sepa- 
rated each by a factor of eight in the particle mass (see 
Table 2). As can be seen, convergence is achieved at 
about hg- The arrows in Figs. 2 and 3 indicate hs for 
the corresponding simulations. In Fig. 3 , the solid-line 
arrow is for the highest-resolution simulation (1024^ ef- 
fective number of particles) , while the dashed line arrow 
is for the 512^ simulation; the NFW fit is shown only for 
these two cases. 

The inner density profiles of all halos simulated here 
are systematically steeper than the corresponding fitted 
NFW law. The slopes of the profiles at r w 1% the virial 
radii, R^,, span a range from —1.4 to —1.6. For compar- 
ison, the slope of the density profile of a typical LCDM 
halo of 2 X IO^'^H-^Mq at O.OIR^ is « -1.2 (a NFW 
profile was used with the corresponding concentration 
given by Bullock et al. 2001, and re-scaled to as = 0.8). 
At radii smaller than O.OIR^,, the slopes tend to become 
shallower but the halos are still denser and slopes steeper 
than the corresponding NFW fit up to the resolution lim- 
its (« hg, see Figs. 2 and 3). 

The overall density profile shapes of our halos are also 
somehow difierent from the NFW function. For radii 
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Fig. 2. — Upper panels: Density profiles of lialos formed at the scale of damping (truncation). From the left, the first two panels present 
the profiles (thick lines) of two halos each one; the lower profiles were shifted by —1 in the log for presentation purpouses. Thin short-long 
dashed lines show the corresponding NFW fits. The third panel shows halo K simulated with four different resolutions. The NFW fits are 
showed only for the two highest resolutions. The innermost radius of each profile corresponds to « 1.3 /14 and arrows in each panel indicate 
the radius /ig (see text). Lower panels: Fractional residuals of the NFW fit. We observe that all our halos are denser and steeper than the 
NFW model at least until ftg. 



larger than 0.02R.U, the profiles tend to be in general 
slightly less curved than the NFW model. This is why the 
residuals shown in Figs. 2 and 3 indicate a systematical 
defect at intermedium radii and then an excess at the 
outer radii. The outer slopes are > —3. 

In summary, the density profiles of the halos simulated 
here, with masses close or below the filtering mass, have 
a shape slightly flatter than the NFW law for radii > 
0.02R^ and slopes significantly steeper at radii < 0.02R„. 
Nevertheless, each profile is different. The profile of halo 
E512o.o has minimal deviations from the NFW function, 
while the profile of halo D512o.o significantly deviates 
from this function. 

Table 1 details the main properties of the halos studied 
here. In columns 6 to 14 are reported respectively the 
virial mass, M^,, the maximum circular velocity, Vmax, 
the Ci/5 and NFW concentration parameters, the aver- 
age density within 1% the virial radius, and three core 
radii estimated by different criteria (the latter quantities 
apply only to halos simulated with Urms, see §§4.2). The 
NFW and C1/5 concentrations are defined respectively as 
the ratios between and the NFW scale radius, and be- 
tween R„ and the radius where 1 / 5th of M^, is contained 
(Avila- Reese et al. 1999). As found in previous results, 
halos of masses below the truncation mass in the power 
spectrum tend to be less concentrated than ACDM ha- 
los of similar masses (Avila- Reese et al. 2001; Bode et 
al. 2001; Knebe et al. 2002). We have simulated some 
ACDM (erg = 0.8) halos of masses « 2 x lO^^ /j-iMq 
and measured NFW and Ci/5 concentrations around 8- 
15 and 6.0-11.0, respectively, to be compared with the 
values given in Table 1 for the Urms=0 cases. 

4.2. The inner structure of halos simulated with rms 
velocity added to the particles 

We have rerun the halos presented in the previous sub- 
section but now introducing a random velocity compo- 
nent, Wrms- As explained in §2, the particle DF used 




Fig. 3. — Upper panel: Density profile for halo D. This halo 
was simulated with four different resolutions, separated each by a 
factor of eight in particle mass (see Table 2). Shown with a short- 
long-dashed curve is the NFW fit to the density profile for the most 
well resolved case (solid line) while arrows indicate the position of 
the radii hg for the two highest resolution simulations. 



corresponds to a Fermi-Dirac function. Regarding the 
v^^s amplitude, we use two values: ■i;rms(0) = 0.1 and 
0.3 km/s. These values are for WDM particles of thermal 
origin and for non-thermal sterile neutrinos, respectively, 
in both cases with mw=0.5 keV. Our goal is to explore 
whether the inner structure of the halos becomes affected 
significantly or not by adding Wims- 

We first present results for the case Wrms(O) = 0.1 km/s, 
and then explore how the inner halo structures change 
when Wrms(O) is increased from 0.1 to 0.3 km/s, a more 
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Fig. 4. — Upper panels: Comparing the density profiles of simulated halos without (solid lines) and with (dot— dashed lines) adding 
''rms to the particles. In the left panel, halos B and B„ have been shifted by —1 in the log. The innermost plotted radii are at 1.3 /i4 and 
the arrows indicate the corresponding hg radii. For halo D, the profile of the 1024'^ particles simulation (nrms=0) is also plotted (dotted 
line). Medium panels: The corresponding 3-D velocity dispersion profiles of the halos shown in the upper panels. Halos B and B^ have 
been shifted by —0.5 in the log. Lower panels: Coarse-grained phase-space density profiles corresponding to the halos shown in upper and 
lower panels. 



appropiate value for the 0.5 keV sterile neutrino model 
used to generate the initial power spectrum of the sim- 
ulations (see §2). For halo D, we have simulated with 
Vims > the 512^ case. Unfortunately, due to limita- 
tions in our computational resources, it was not possible 
to run the simulation with 1024^ particles for v^ms > 0. 

In Fig. 4 we present spherically-averaged profiles of 
different properties for all the simulated halos without 
and with Vims added. In the first and third columns 
two halos (A and B, and Db and D, respectively) are 
presented but the latter ones are down shifted for clar- 
ity. The upper panels of Fig. 4 show the density pro- 



files of our halos without (black solid lines), and with 
Wrms(O) = 0.1 km/s (red dot-dashed hue) and 0.3 
km/s (magenta dashed line). In the medium panels are 
plotted the corresponding three-dimensional velocity dis- 
persion profiles, cr3D(''), and the lower panels show the 
corresponding coarse-grained phase-space density pro- 
files, Q{r) = p{r)/a^j^. As in Fig. 2, the arrows indicate 
the strong resolution limit radius /ig of the simulations. 
For halo D (third column) , the simulations with v^ms > 
were carried out with 512^ particles and using two dif- 
ferent random seeds for the particle DF calculation (see 
below). 
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To some degree, all the simulated halos have been af- 
fected in their inner regions by the injection of initial ran- 
dom velocities to their particles. For the less resolved ha- 
los A, B (left panels) , the deviations at ~ hg of the inner 
density profiles from the profiles obtained in the simula- 
tions with i'rms=0 are yet marginal, but in the expected 
direction. For the halos C, D, Db, and E simulated with 
512^ particles, the deviations down to /ig are significant: 
the density profiles systematically flatten with respect to 
the corresponding fims=0 cases. The radii at which the 
density profiles of halos with Vrms start to deviate (flat- 
ten) from the ones without Urms, are ^ 0.010 — O.OlSRi,, 
well above the resolution limit of the simulations. 

In a second series of experiments, we have resimulated 
halos C, Db and E (512'^ particles) with tir,„s three times 
larger, i.e. frms(O) = 0.3 km/s. The profiles correspond- 
ing to these simulations are plotted in Fig. 4 with dashed 
(magenta) lines. The flattening of the inner density pro- 
files is clearly more pronounced than in the simulations 
with Wrms(O) = 0.1 km/s. 

The innermost density proflles actually vary from halo 
to halo. Again, halo E is the less affected not only by 
the damping in the power spectrum, but also by the 
injection of Wrms, and halo D is the most affected by 
both effects (lower dot-dashed curve in the correspond- 
ing panel). The latter actually shows a "true" flat core 
already at h^. Since the core of this halo is too different 
with respect to the other ones, we decided to explore if 
large difference can be explained by a rare fluctuation 
in the random procedure of particle velocity assignment. 
Thus, the same halo D512o.i was resimulated with a 
different seed in the random number generator used to 
draw the particle velocities. The upper curves in the 
third column panels of Fig. 4 correspond to the proflles 
for this halo, called D512o.ib. The inner density proflle 
of this halo is not too different from the proflles of the 
other halos, though it remains as the flatest one among 
all the simulated halos with Wrms(O) = 0.1 km/s. 

In Fig. 5 we attempt to flt different functions to the 
density proflles of halos C, D, and E (512^ particles) with 
'^rms(O) — 0.1 and 0.3 km/s. A general function to de- 
scribe density proflles of cosmic objects was proposed by 
Zhao (1996): 



' (r/ro)'/[l + (r/ro)"](/3-^)/" ' ' 

The NEW proflle corresponds to {a,(3,j) = (1,3,1). 
This function does not provide a good description for 
the proflles of our halos with Vi-ms > 0, in particular in 
the inner regions. We have fltted the halo proflles to 
the NEW in order to obtain an estimate of the Cnfw 
concentration given in Table 1. 

Strigari et al. (2006) suggested a "cored" density 
proflle in order to derive constraints on the size of a 
possible shallow core in the halo of the Fornax dwarf 
spheroidal galaxy (Goerdt et al. 2006; Sanchez-Salcedo 
et al. 2006). This proflle is described by eq. (5) with 
{a,(3,j) — (1.5,3,0), and they deflne the core radius, 
ic as the radius where the inner slope, g, reaches the 
value of -0.1. Thus, r^ = ro/{-3/g - 1)^/" « O.lro. 
The dot-dashed curves in Fig. 5 show the best flts us- 
ing the S2006 proflle. As can be seen, this proflle does 
not describe well the density proflles of the WDM halos 
in the simulations. Note that a characterizes the sharp- 




FlG. 5. — Density profiles of simulated (512'^ particles) halos with 
added Drms and fitted to different model profiles. Long-dashed 
line: {a,/3,7) = (0.7,3,0); dot-dashed line: the profile proposed 
in S2006, (o,/3,7) = (1.5,3,0); dotted line: the same S2006 profile 
but fitted only to the inner 30 kpc. 



ness of the change in logarithmic slope. As already seen 
in Figs. 2 and 4, the proflles of our halos tend to be 
less curved than the usual NEW proflle. Therefore, val- 
ues of a smaller than 1 should be used instead of larger 
than 1. We have obtained a reasonable description of our 
WDM proflles with (a,/3,7) = (0.7,3,0). The (magenta) 
dashed curves in Fig. 5 are the best flts with these pro- 
files. The core radius, as defined above, is in this case 
Tc ~ 0.0064ro. Finally, we have also tried flts to the 
S2006 function but taking into account only the central 
halo regions, up to ~ 30/i~^kpc. The flts are shown in 
Fig. 5 with (green) dotted curves. The flt is specially 
good for the halo D512o.i and those with Wrnis(0)=0.3 
km/s. Columns 12 to 14 in Table 1 report the values of 
Tc obtained with the three different flts. 

As to the 3-D velocity dispersion proflles of the simu- 
lated halos, they do not differ signiflcantly between the 
cases with and without Wrms, the exception being halo D. 
In the innermost regions, i73d('') is similar or higher for 
halos with Vrms than for those without fims- The larger 
differences are for halo D, which after introducing ran- 
dom velocities to the particles produce a relatively hot 
core. 

Finally, from the measured density and dispersion ve- 
locity proflles, we calculate the coarse-grained phase- 
space density proflles, Q{r). As seen in the lower pan- 
els of Fig. 4, excepting the innermost regions, the Q{r) 
proflles are well described by a power law Q cx r^" with 
a ~ —1.9 close to that obtained for CDM halos by Taylor 
& Navarro (2001). For the inner regions, the Q{r) pro- 
flle of the halos simulated without Vrms tends to steepen. 
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specially in halos D512o.o and C512o.o- The opposite 
happens for the halos simulated with adding Urms, the 
inner Q{r) profile tends to be flatter as Vims is larger. 

5. DISCUSSION 

5.1. Robustness of the results 

The halos studied here are on one hand among the first 
virialized structures to form in our simulations, and on 
the other their assembling process started relatively late 
in the universe, between z « 0.6 and 1 (for discussions 
about the mass assembling process of halos with masses 
close to the damping scale in the power spectrum see e.g., 
Moore et al. 2001; Avila-Reese et al. 2001; Bode et al. 
2001; Knebe et al. 2002; Busha et al. 2006). Because of 
the late collapse, one might argue that the halos studied 
here arc not relaxed and therefore is not suprising to have 
the deviations from the NFW profile, as reported in §4.1 
for the experiments with i)i.nis=0. We have followed the 
evolution of some of our halos for more than a Hubble 
time (scale factor a > 1), finding a negligible evolution 
in the density profiles since a ~ I. For example, halo 
D512o.o was run until a = 1.4 (18.4 Gyr). The density 
profile of the halo at this epoch is practically the same 
as at a = 1 (13.7 Gyr, see Fig. 4). As discussed in 
previos studies (see the references above), the collapse of 
halos of scales close to the damping scale seems to be 
quasi-monolithic (though highly non spherical). Thus, 
in regions that remain relatively isolated, as in the case 
of the halos selected for our study, halos suffer low mass 
accretion, and their structures remain almost unaltered 
since the initial collapse. 

We checked that the effects on the structure of halos 
reported in §4 are systematic by running for a WDM 
model -not shown in Table 1- the corresponding CDM 
simulation, using the same random phases and changing 
only the initial power spectrum. We found that the den- 
sity profile of the CDM halo is well fitted by the NFW 
function, while the corresponding WDM halo present the 
systematic deviation already seen in Fig. 2. Figure 6 
compares the density and circular velocities profiles for 
the halo in question in its two versions, CDM and WDM 
without random velocities. Since the WDM halo ends up 
at ^; = with a slightly lower mass than the CDM halo 
we correct the profiles so as to make the comparison at a 
fixed mass (= 3.0 x 10^^ h~^MQ). Notice, in particular, 
that the inner density profile of the WDM halo is indeed 
steeper than their CDM counterpart. 

Concerning the resolution limit in our simulations, 
based on the convergence study carried out for halo D 
(see Fig. 2), we find that a strong limit is /ig = S/ifor, 
but resolution might be still acceptable for radii slightly 
larger than /14 = 4/ifor, a value suggested previously for 
CDM halos in "equilibrium" simulated with the ART 
code (Klypin et al. 2001). With a resolution limit at 
hg, our simulations allow to resolve the inner structure of 
halos down to 1.2/i~-'^kpc for the 512^ runs and down to 
0.61/i~^kpc for the 1024"^ run. These radii correspond 
respectively to ^ 0.5% and 0.25% virial radius in our 
halos. 

Finally, it is important to recall that the masses of the 
the WDM halos analyzed in our simulations, are well 
above from the mass scale affected by discreteness ef- 
fects, like the spurious formation of structures and sub- 



structures due to the initial grid pattern (see §3.1). 

5.2. Do soft cores form in WDM halos? 

Early structure formation studies based on a WDM 

cosmology considered particles originated in thermal 
eqiiilibrium. For this case, both ?;rms and A j depend 
only on the particle mass tow- The smaller mw, the 
larger A / and Vrms ■ Controlled numerical simulations of 
isolated halos showed that in order to produce "observ- 
able" soft cores, the amplitude of Wrms should be several 
times higher than the values corresponding to thermal 
WDM particles of masses mw ^ 1 keV (Avila-Reese et 
al. 2001). Particle masses smaller than ~ 1 keV are not 
allowed by the constraints on satellite galaxy abundances 
as well as by the Ly-a power spectrum alone or combined 
with CMBR and large scale structure data. The Ly-a 
power spectriim is the strongest of the constraints. For 
the non-thermal sterile neutrino, it places a limit on it 
mass at mw > 2 keV (Seljak et al. 2005; Viel et al. 2006; 
Abazajian 2006b). For thermal WDM particles, the ob- 
servational constraints give a limit of mw 0.5 keV 
(Narayanan et al. 2000; Viel et al. 2005; Abazajian 
2006b), while a distinct analysis, using different simula- 
tions provide a stronger limit, mw ^ 2.5 keV (Seljak et 
al. 2006). 

We may estimate the expected flat core radii of 2 x 
10^'^h~^MQ WDM halos for thermal particles in the mass 
range mw = (2 — 0.5) keV by using the approximation 
given in Avila-Reese et al. (2001, their eq. 13). This ap- 
proximation is based on the monolithic collapse of halos 
with non-negligible particle random velocities before the 
collapse. For thermal WDM particles of masses 2 and 0.5 
keV, the zm of a 2 x IO^'^H^^Mq perturbation are « 3.4 
and 1.9, respectively, while the Urms(O) corresponding to 
these masses are 0.1 and 0.015 km/s. Therefore, the ex- 
pected core radii are rc w 30 and 210 pc. 

The value of Wi.,„s for a non thermal sterile neutrino 
of 0.5 keV is approximately three times larger than the 
corresponding to the thermal particle of the same mass. 
Therefore, for this case ic ~ 630 pc. So, the resolutions 
that we may attain in our simulations of WDM halos 
for the mw=0.5 keV sterile neutrino are already close to 
these estimates of the flat core radii. 

Recently, alternative particle models, like super- 
WIMPS, were proposed. The dark particles in these 
models may acquire random velocities non-thermally, 
for example, through the decay process of NLSP par- 
ticles (e.g., charged Sleptons into gravitinos, Feng, Ra- 
jaraman & Takayama 2003, and see for more references 
Steffen 2006). In these does not depend di- 

rectly on the damping scale of the linear power spectrum. 
However, an extra parameter is introduced, the decaying 
epoch. Strategies to constraint this parameter using as- 
trophysical observations have been proposed (Feng, Ra- 
jaraman & Takayama 2003; Strigari et al. 2007). 

Some cosmological models with super- WIMP particles 
may be in agreement with constraints based on structure 
formation, specially the Ly-a forest, and still allow for 
relatively large velocity dispersions, able to ameloriate 
the potential problems of ACDM at small scales. This 
is the case of the so called Meta-dark matter models 
that consider the late decay of neutralinos into graviti- 
nos. These models preserve a power spectrum similar to 
ACDM models and at the same time set a phase space 
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limit in the innermost structure of dark halos due to the 
injection of random velocities to the particles (Strigari 
et al. 2007; Kaplinghat 2005). However, as mentioned 
in the Introduction, it was still an open question of how 
efficient the introduction of Wmis is for producing sig- 
nificant effects on the inner structure of simulated dark 
halos. One of the goals of the present paper was just to 
expore this question by means of numerical simulations 
able to resolve the WDM halos down to ~ 0.005 Ri,. 

The results presented in §4.2 show definitively that 
the addition of Urms flattens the inner density profiles 
of WDM halos, and more as Vrms is higher. Differ- 
ences in density profiles of halos simulated with Wrms > 
and those with Vi-ms = 0, start to be evident at radii 
of 2-3 /i-ikpc(- O.OlORi, - 0.015R„) for our w 2 - 4 x 
10^^/i~^Mq halos. The halo average inner density mea- 
sured at O.OlR^u, pi% decreases on average by factors 1.5 
and 2.5 for the models with Wrms(O) = 0.1 and 0.3 km/s, 
respectively (column 11 in Table 1). Nevertheless, we 
notice that the concentration of the halos seems not to 
change significantly by the introduction of Urms- 

In general, the NFW function does not describe well 
the inner density profiles of our simulated halos; their 
profiles are much shallower than 7 = — 1 as can be seen 
in Fig. 5. However, is the size of the random velocity 
effect as theoretically expected? The theoretical predic- 
tions are based on the existence of an upper limit in the 
fine-grained phase-space density due to the coUisionless 
nature and finite relict velocity dispersion of particles. 
This upper limit, Qo.max, implies that the halo density 
profile must saturate and form a constant-density core 
(Hogan & Dalcanton 2000). 

For the random velocities used in this study (wrms(O) = 
0.1 and 0.3 km/s, corresponding respectively to thermal 
and non-thermal 0.5 keV sterile neutrinos), Qo.max = 
3 X 10-5 and 1.1 x 10"^ MQpc-3/(km/s)3. According 
to Hogan & Dalcanton (2000), the core radius produced 



by the phase-space packing scales as m 



-1/2 



X 



where 



f c,oo is the asymptotic circular velocity for the assumed 
non-singular isothermal sphere (their eq. 18). This im- 
plies that more massive halos have smaller, more tightly 
bound cores. Applying the same equation for mx — 0.5 
keV and Vcoo = 200 km/s, the core radius for the ther- 
mal particle (Qo,max = 3 x 10^*^ M0pc-^/(km/s)^) is 
85 pc, while for the sterile neutrino (Q = 1.1 x 10^^ 
MQpc-3/(km/s)3) is 450 pc. 

In any case, our results can only marginally test these 
estimates. The resolution limit in our simulations with 
non-zero Urms is in between w 1 and 1.7 kpc. If the inner 
density profile response to the introduction of ?;inis is 
gradual, one expect to see yet some effects at these radii, 
and this happen to be the case. 

A way to attempt to infer (extrapolate) the sizes of 
possible flat cores in the simulated halos is by fitting the 
measured density profiles to an analytical function that 
implies a fiat core. Results of these fits were presented 
in §4.2 using the Zhao (1996) profile with (a,/3, 7) = 
(0.7, 3, 0) as well as the one suggested by S2006. The lat- 
ter function gives a poor description of the overall mea- 
sured density profiles, which tend to be significantly less 
curved than the analytical model (see Fig. 5). The ob- 
tained (overestimated) values for {— O.lrg, see §4.2) 
are reported in column 13 of Table 3. When the S2006 
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Fig. 6. — Upper panel: Comparison of the density profile of a 
WDM Milky Way-like halo, without the random velocity, (squares) 
with the density profile of its counterpart CDM, this latter gener- 
ated using the same random phases (crosses). Lower panel: Same 
as upper panel but for circular velocity. 



function is fitted to only the inner 30 kpc, the fits im- 
prove and the estimated core radii become smaller by 
roughly a factor of two (column 14). However, even for 
this case the core radii seems to be upper limits, with 
the exception of halo D512o.i. 

We have found that the WDM profiles are better de- 
scribed by the Zhao function with (a,/3, 7) = (0.7,3,0). 
The best fits to the measured profiles give extrapolated 
core radii rc « 5 — 8 times smaller than the S2006 profile 
fitted to only the inner 30 kpc. Should we have the suffi- 
cient resolution to resolve the flat cores, their radii would 
lie in between Vc and ?'f 30- The only halo for which the 
flat core is patent at our resolution limit is D512o.i; a 
visual inspection shows that the core radius is close to 1 
h~^kpc. 

Our results show that the different estimates of the 
core radius increase by less than a factor of 1.5 from 
the simulations with '(;„„s(0)=0.1 km/s to the ones with 
Wims(0)=0.3 km/s. The amount of this increase is less 
than the one we would predict using the monolithic col- 
lapse approximation of Avila-Reese et al. (2001); ac- 
cording to this approximation, the core radius of halos 
formed at the same time depends linearly on the injected 
Uims at the maximum expansion of the perturbation, 
ic oc Wrms,ZM = ^'rms,o(l + zm)- We have estimated the 
predicted values of Tc for our halos by using this approx- 
imation. From the simulations, we find that the redshifts 
of maximum expansion of halos C, D, and E are roughly 
zm = 1-6, 1.3 and 1.8, respectively; these redshifts are 
practically the same for the different values of Urms- The 
calculated Tc for Wrms(O) = 0.1 (0.3) km/s are then 221 
(663), 240 (720), and 209 (627) pc, respectively. Thus, in 
general these predictions give core radii just in between 
Tc and rfgQ (see Table 1), though the dependence on 
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frms is much more pronounced than for the (extrapo- 
lated) core radii estimated from the fits to our halos. 

5.3. How is the structure of halos at the damping 

scale? 

The simulations carried out in this paper allowed us 
also to explore the; structure of dark matter halos formed 
from perturbations at the scale of damping of the lin- 
ear power spectrum. The cutoff in the power spec- 
trum used here corresponds to a relatively large mass, 
Mf = 2.6 X lO^^/j-i^-ig Therefore, the formation of 
the first structures in this model, namely those struc- 
tures with masses close to Mf, happens relatively late. 
We speculate that the formation process of the trunca- 
tion halos is generic. If this is true, then the structure 
of the late-formed truncation halos simulated here with 
several millions of particles (wrms=0) should be similar to 
the structure of early-formed truncation (micro)halos in 
models with much smaller filtering masses than the one 
used here, for example in the CDM models. If this is the 
case, then our results may enrich the discussion about the 
formation and structure of the first microhalos in CDM 
models (Earth-mass scales). 

We have found a clear systematic trend in the den- 
sity profiles of the simulated truncation halos: they 
are significantly steeper than r~^ in the inner regions, 
r < 0.02R,,, - 0.03R„, and he below the best NFW fits in 
the intermediate region. CDM halos with an inner slope 
steeper than r~^ have been reported in other contexts: 
recently merged group- and cluster-size halos (Knebe et 
al. 2002; Tasitsiomi et al. 2004) or microhalos formed at 
the scale of CDM power spectriim damping (Diemand 
et al. 2005). It has been argued that the recent ma- 
jor merger is the one to blame for the steepening of the 
density profile while subsequent secondary infall modifies 
the external region. For the halos at the damping scale, 
rather than a major merger, the dynamical situation cor- 
responds to a fast (cuasi-monolithic) collapse. However, 
in both cases, the process is dynamically violent. We 
have checked that the obtained density profiles do not 
correspond to a transient configuration. As mentioned 
in §5.1, for halo D the profile remains almost unchanged 
until a = 1.4. 

Extrapolating our results to the damping scale of 
CDM, we can speculate that CDM microhalos may be 
significantly steeper than the NFW profile. This im- 
plies that the possible contribution of surviving micro- 
halos to the 7— ray fiux originated by neutralino anni- 
hilation, might be comparable to the central fiux from 
host halos (Diemand et al. 2006). Another implication 
of our results could be related to the buildup of the inner 
density profile of dark halos in general. Dehnen (2005) 
and Kazantzidis et al. (2006) argued that the assembly of 
halo inner density profiles happens very early in the his- 
tory of the universe and specifically that the inner slope 
of the cuspier progenitor survives up to the final halo. If 
the density profiles in our simulated halos are represen- 
tative of objects formed at the damping of power scales 
in general, then, according to the mentioned studies, the 
central slope of present day dark matter halos should be 
much steeper than the NFW profile. 

In order to verify our results, a more systematic study 
halo structure at the scale of damping is required, ex- 
ploring any possible dependence with the shape of the 



cutoff and the power spectrum slope at the scale of damp- 
ing. Currently the only studies discussing a similar situ- 
ation report different results for the scale of galaxy clus- 
ters(Moore et al. 1999b) and the microhalos (Diemand 
et al. 2005). It is unclear if this suggests a slope depen- 
dence on the profile of the smallest dark matter halos 
with the power spectrum slope. 

6. CONCLUSIONS 

We have studied by means of cosmological N~body 
simulations the structure of dark halos formed in the con- 
text of a WDM model corresponding to a non- thermal 
sterile neutrino particle of mass mw=0.5 keV. The first 
series of simulations did not include the injection of a 
random velocity component, Vrms, to the particles and 
were aimed at exploring the structure of the halos formed 
from perturbations at the damping scale in the linear 
power spectrum (Mf = 2.6x10^^ h~^MQ for the concrete 
WDM model studied here). The second series of simula- 
tions included a Wrms component of a (i) a mw=0.5 keV 
thermal neutrino (vrmsiz = 0) = 0.1 km/s), and one that 
roughly corresponds (ii) to a mw=0.5 keV non-thermal 
sterile neutrino {vrmsiz = 0) « 0.3 km/s). This latter 
model was used to generate the initial power spectrum. 
These simulations were aimed at exploring the effect that 
a random velocity component has on the inner structure 
of halos; in particular, at dilucidating whether constant 
density cores are produced or not in the simulations. The 
results of our study lead us to the following two main 
conclusions: 

• The structure of halos formed from perturbations of 

scales close to Mj, and resolved with up to more than 
16 millions of particles (with Urms=0), is peculiar: the 
inner density profile (< 0.02R„) is systematically steeper 
than the best corresponding NFW fit (and the respec- 
tive CDM counterpart), and the overall density profile 
(> 0.02R„) tends to be less curved than the best NFW 
fit; the outer profile slope is never steeper than —3. Ac- 
cording to our tests, these differences with respect to the 
structure of halos assembled hierarchically, can hardly 
be attributed to a peculiar dynamical state of the halos 
simulated here. 

• The effect of adding Wrms to the particles pro- 
duces a significant flattening of the inner density profile 
(r < 2-3/i-ikpc corresponding to ~ O.OIOR^ -0.015R„) 
of the simulated halos. The different estimated (extrap- 
olated) sizes of the nearly constant density cores are of 
the order of the theoretical predictions, which give values 
below our resolution limit. For the halo masses simulated 
here, M^ « (2 - 4) x 10^^ /i-^M©, the fiat core radn es- 
timated from different fittings are between ~ 0.1 — 0.8 
h~^]ipc. An increase in Wrms(O) from 0.1 to 0.3 km/s pro- 
duces an increase in the extrapolated core radii of a fac- 
tor 1.5 or less. For one of our simulations (halo D512o.i), 
the presence of a nearly constant-density core, of radius 
w 1 /i~^kpc, is already revealed at the resolution limit; 
the same halo simulated with a different random velocity 
seed is less flattened. 

Although the simulations presented here refer to a 
concrete WDM model, they can be interpreted within 
a wide range of contexts. The density profile of dark 
halos with masses close to the truncation mass in the 
linear power spectrum is systematically different from 
the NFW profile; in particular, the inner regions tend 
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to be steeper. These eould have important impUcations 
in the context of CDM models if a significant fraction 
of microhalos formed at the free-streaming CDM scales 
lO~^/i~^M0) have survived until the present epoch. 
In this case, the predicted 7— ray flux from the neu- 
tralino annihilation in the center of these cuspy micro- 
halos might be comparable to the central flux from host 
halos. On the other hand, the fact that the microhalos 
are so cuspy, could have some interesting implications in 
the building up of the next hierarchies of the halo as- 
sembling, as well as in the inner structure of the larger 
halos. 

Regarding the effects of random velocity injection to 
the particles, our results show evidence of significant in- 
ner flattening of the halo density profile at our resolution 
radii. These resolutions are not enough to test directly 
the predicted core radii by phase space constraints (e.g., 
Hogan & Dalcanton 2000) or by dynamical models of 
gravitational collapse with initial random velocities (e.g., 
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